Direct X-ray Constraints on Sterile Neutrino Warm Dark Matter 
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Warm dark matter (WDM) might more easily account for small scale clustering measurements 
than the heavier particles typically invoked in A cold dark matter (ACDM) cosmologies. In this 
paper, we consider a AWDM cosmology in which sterile neutrinos v B , with a mass m s of roughly 
1-100 keV, are the dark matter. We use the diffuse X-ray spectrum (total minus resolved point 
source emission) of the Andromeda galaxy to constrain the rate of sterile neutrino radiative decay: 
v s — > Ve,fi.,T + 7- Our findings demand that m s < 3.5 keV (95% C.L.) which is a significant 
improvement over the previous (95% C.L.) limits inferred from the X-ray emission of nearby clusters, 
m s < 8.2 keV (Virgo A) and m s < 6.3 keV (Virgo A + Coma). 

PACS numbers: 95.35.+d, 13.35.Hb, 14.60.St, 14.60.Pq 



I. INTRODUCTION 



Standard ACDM models of galaxy formation predict 
more small scale structure than is observed, e.g., ex- 
cess galactic satellites ("the missing satellite problem" 
0, H IE 3), high galactic central densities ("the cen- 
tral density problem" |E IE II IE H)> etc - Because 
structure formation is suppressed on scales smaller than 
Xfs — 0.25 Mpc (mW keV)~ 4 ^ 3 for dark matter par- 
ticles of mass mvM HE] > models of Warm Dark Matter 
(WDM) with ~ keV masses, such as sterile neutrinos 
[Til Il2l Il3| . may lead to better agreement with obser- 
vations than standard (GeV - TeV) CDM candidates, 
like the neutralino. On the other hand, as much as the 
suppression of small scale structure by WDM may al- 
leviate problems at low redshifts, the same suppression 
will also delay the onset of reionization at higher red- 
shifts HI [13, possibly to an extent that is difficult or 
impossible to reconcile with the WMAP3 results (lt| . 
Constraints on the level of small scale clustering from 
measurements of the CMB H EE 13 HE El El , the 
3-D ga l axy power spectrum |23| , and the Lyman-a forest 

EE EE El EE EE HI also P reclude drastic small scale 
suppression. In the context of the sterile neutrino (v s ) 
WDM model presented in [H |H |H E| , these data re- 
quire m s > 1.7 keV (95% C.L.) according to the analysis 
presented in EE]- More recent work by Seljak et al. |34j ] 
suggests that this lower bound can be improved by al- 
most an order of magnitude to m s > 14 keV (95% C.L.). 
Because indirect limits imposed by small scale clustering 
data require the interpretation of simulations at their res- 
olution limit and because of the striking improvement of 
this constraint over previous work, this result has already 
generated some debate [33, 36] . Even more recently, how- 
ever, Viel et al. |35| reported a very similar indirect limit 
of m s > 10 keV (95% C.L.), essentially confirming the 



results of Ref. [3J. 

Fortunately, it is also possible to directly constrain m s 
based on the radiative decay of sterile neutrinos to X-rays 



of energy E 7 ^ s — m s /2 via v s 



-7. 



This test is di- 
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rect in the sense that it probes the signature of individual 
particle decays, whereas the cosmological tests are indi- 
rect in the sense that they probe only the macroscopic 
clustering signatures of the dark matter. In order to pre- 
vent a potentially viable dark matter candidate from be- 
ing dismissed prematurely, it is important to separately 
improve both constraints. This is especially true in the 
case of sterile neutrinos which, apart from their possible 
role as the dark matter, have interesting implications for 
several important physical processes including the pro- 
duction of the baryon [H EE EE EE and Jepton H3 
asymmetries, Big Bang Nucleosynthesis Ell |42|, r eion- 
ization [i l El [45I [il neutrino oscillations [TlMcl l49 | . 
pulsar kicks lHaMlHl EE El, etc. 

Several direct limits on m s have already been re- 
ported. The absence of anomalous line features in XMM- 
Newton (0.5-12 keV) EE El and HEAO-I (3-60 keV) 
El E| measurements of the Cosmic X-ray Background 
(CXB) allowed Boyarsky et al. EE to set a (95% C.L.) 
upper limit of m s < 9.3 keV (assuming, as with all the 
mass limits quoted below, that Qdm = = 0.24). In 
Ref. E3, Abazajian, Fuller, and Tucker (hereafter AFT) 
suggested that the best constraints could be achieved by 
examining individual objects, e.g., galaxies and/or clus- 
ters, rather than the CXB. Based on XMM-Newton ob- 
servations of Virgo A (M87) , the dominant galaxy in the 
northern part of the Virgo cluster EE], Abazajian EE 
arrived at an improved upper limit of m s < 8.2 keV 
(95% C.L.). Boyarsky et al. E3 have also used XMM- 
Newton observations of Virgo A EE] and the Coma 
|62ll63| cluster to explore possible constraints on m s , but 
they do not provide a definite numerical limit. Abaza- 
jian and Koushiappas EE] estimate an upper bound of 
m s < 6.3 keV (95% C.L.) based on the results in Ref. [61] . 

An important general feature of all sterile neutrino 
mass constraints is that m s is very weakly dependent 
on the flux from which it is inferred: m s oc EE E| 
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(Eqn. H3 below) . Consequently, even a significantly im- 
proved constraint on <E> X s leads to only a modest im- 
provement in the mass limit. This dependence puts a 
premium on carefully selecting an object with the best 
ratio of dark matter decay signal to astrophysical back- 
ground. 

In this paper, we consider the limits imposed by XMM- 
Newton observations of the Andromeda galaxy (M31) 
[64j . Although these observations probe significantly less 
dark matter than observations of clusters, the predicted 
sterile neutrino decay signal from Andromeda is com- 
parable to that of Virgo A due to Andromeda's close 
proximity (Dm31 — 0.78 ±0.02 Mpc |66j. as compared to 
Z?M87 — 15-8 ± 0.8 Mpc 67]). Moreover, by analyzing a 
single galaxy rather than a cluster, we avoid one substan- 
tial source of astrophysical background: hot, intra-cluster 
gas. Shirey et al. 64] eliminate a second major astro- 
physical background by identifying and removing the flux 
from resolved X-ray point sources in Andromeda. The 
only remaining astrophysical contributions to the diffuse 
spectrum of Andromeda are unresolved point sources and 
(a modest amount of) hot gas. The analysis of Shirey et 
al. [H3] suggests that hot gas dominates the diffuse emis- 
sion at low energies but falls off rapidly at E~ ^ 0.8 keV 
until the unresolved point sources begin to dominate at 
E 1 > 2 keV. This fact underscores the importance of us- 
ing high resolution spectra to constain m s , without which 
we could not benefit from the rapidly decreasing astro- 
physical emission in the energy range of interest {E 1 ^ 
0.85 keV) and would, instead, be forced to compare v a 
decay peaks to much larger broad-band (e.g., soft band: 
0.5-2 keV) fluxes. For these reasons, the diffuse spectrum 
of Andromeda is an excellent data set for constraining 
the decay of sterile neutrinos and, in particular, demands 
m s < 3.5 keV (95% C.L.), as we show below. 

In Sec. m we provide a brief overview of the ster- 
ile neutrino WDM model [3 EH. In Sec. Ell we de- 
scribe the XMM-Newton observations of Andromeda [64| 
and Virgo A [gjj, an d other relevant physical properties 
of these systems (see Table I). In Sec. IIVI we explain 
our analysis and compare the potential for direct detec- 
tion of keV dark matter decay in nearby galaxies and 
clusters. In Sec. we present an updated sterile neu- 
trino exclusion plot, including our mass limit based on 
the diffuse spectrum of Andromeda as well as the previ- 
ous direct and indirect limits. We conclude in Sec. IVII 
Throughout the paper, we assume a flat cosmology with 
^baryon = 0.04, fi WDM = « s = 0.24, ft A = 0.72, and 
h = iJo/100 km s" 1 Mpc" 1 = 0.72. 



light, active neutrinos v a (i.e., v c ^ tT ) at a rate of 
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where 9 is the vacuum mixing angle for an effective 2x2 
unitary transformation from v s <-> v a ; since we consider 
only very small mixing angles, v s is nearly a pure mass 
eigenstate 0, E3 • The radiative decay in which we are 
interested, v s — > i/n + 7, is suppressed by a factor of 
27a/8n ~ 1 / 128 [68l relative to and occurs at a 
rate of [H EH 



T s ~ 6.8 x 10 
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While the decay rate is very slow, a large collection of 
sterile neutrinos will produce a detectable X-ray signal. 
In particular, the X-ray luminosity resulting from the 
decay of N s — (Mdm/w s ) sterile neutrinos in a dark 
matter halo of mass Mdm is given by 
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The corresponding line flux at _B 7 S — m s /2 is 



5.1 x 10 18 erg cm 



A/dm 
10 n M Q 
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If we assume a QCD phase-transition temperature of 
Tqcd — 170 MeV and a lepton asymmetry of L ~ 7710 = 
^baryon/^7 — 10 _ 10 ) |32lj3q] . the sterile neutrino density- 
production relationship |32] (updating [3l|]) is 



= 3.27 keV 



sin 2 2(9 

10- 8 



-0.615 



0.24 



0.5 



(5) 



Although the values of Tqcd and L are uncertain, we 
adopt Eqn. for definiteness and for ease of compar- 
ison to the literature. Our limit on m s could easily be 
re-evaluated for a different density-production relation- 
ship; for example, note the three lines for models with 
large L in Fig. |21 below. By combining Eqns. (@J and 
we arrive at an expression for the line flux that is 
independent of the mixing angle: 



II. THE STERILE NEUTRINO WDM MODEL 



In the model presented in Refs. |12|, |31j, |32j, |33| , sterile 
neutrinos v s with m s < m e predominantly decay to three 



$x,s(^ S ) 



3.5 x 10~ 15 erg citTV 1 f ) 
6 V M P C / 
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3.374 
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In the remaining sections we describe the data and meth- 
ods we use to constrain m s in the context of this model. 
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III. PROPERTIES OF ANDROMEDA AND 
VIRGO A 

Comparison of the X-ray emission from Andromeda 
and Virgo A is facilitated by the fact that both sys- 
tems were observed with the same instrument: XMM- 
Newton [60l l6^ | . Even more important is the fact that 
XMM is sensitive to the full (95% C.L) range of presently 
allowed v s radiative decay energies (0.85 keV < E ltS — 
m s /2 < 3.15 keV) j3||. In addition to their X-ray prop- 
erties, this section also provides the distances and dark 
matter mass estimates we adopt for our calculations of 
the sterile neutrino decay fluxes from each object. 



A. X-ray Data 

Shirey et al. [(54| observed Andromeda with XMM for 
34.8 ks out to a radius of 15 arcminutes (~ 3.4 kpc at 
£>M3i - 0.78 ± 0.02 Mpc @). The (0.5-12 keV) diffuse 
spectrum we utilize was extracted from within a radius 
of 5' (~ 1.1 kpc) from the center of the galaxy. Shirey 
et al. 0| produce the diffuse spectrum by identifying 
and removing the flux from resolved X-ray point sources 
in Andromeda; we make no attempt to model or sub- 
tract the remaining astrophysical contributions from un- 
resolved point sources and hot gas. We note that the 
Shirey et al. data are consistent with the (0.5—7 keV) 
diffuse spectra presented in a joint Chandra— XMM study 
of Andromeda by Takahashi et al. [65| . 

Bohringer et al. |(3(j used XMM to measure the flux of 
and around Virgo A (M87) out to a radius of 12' (~ 55 
kpc at £>M87 — 15.8 ± 0.8 Mpc (63) over a (usable) ex- 
posure time of 25.9 ks. Abazajian et al. 0,H3 consider 
the (0.5—8 keV) flux from within a radius of 8.5' (~ 39 
kpc). These data are summarized in Table I. 



B. Dark Matter Masses 

To estimate the dark matter mass of Andromeda 

and Virgo A that is enclosed within the XMM field of 
view (fov), we integrate the (r -2 — weighted) dark matter 
density of each halo pdmC^* - D\) over a truncated cone 
of radius i?f ov and length 2R vlr : 



Sfov — 



PDM(\r- D\)dV [ov 



(7) 



We define = Z? 2 Ef ov . The distance to the center of 

each object is D, i? v ; r is the virial radius, 



Rfov = iov r ~ 0.3 kpc I -rjy- 



?fov 



Mpc 



(8) 



is the radial extent of the fov at a distance r, and r varies 
between D — i? v ; r , at the near "edge" of each dark matter 
halo, to D + i? v i r at the far "edge" . 



Based on rotation curve data, Klypin, Zhao, and 
Somerville |70j estimated the dark matter mass distri- 
bution of Andromeda, pdm,M31- When we integrate 
PDM,M3i over Vf ov ,M3i (0fov = 5' 0; i? vir ~ 300 kpc 
|70| ) , we find that the region from which the diffuse emis- 
sion spectrum of Andromeda was extracted contains 
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JW DM,M31 



(0.13 ± 0.02) x 10 1X M ( 
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As shown in Ref . [7(j , about half of this mass is enclosed 
within a sphere of 1 kpc radius about the center of An- 
dromeda. It is worth re-emphasizing that our mass limit 



scales like J cx (M, 
sitive to uncertainties in Af§^ 



dm) ' 3 ? and is therefore very insen- 



To be conservative, we 
ignore the contribution from the fraction of the Milky 
Way halo within the fov; see Refs. |40j, |7lJ for discussion 
of possible constraints based on Milky Way dark matter 
alone, i.e., "blank sky observations." 

In Ref. [U , AFT integrate an isothermal /3-model 0| 
to estimate the dark matter mass of Virgo A within 
a 17' x 17' rectangular prism. When we integrate the 
same model over a truncated cone of cross-sectional ra- 
dius fov = 8.5' and length 2i? vir ~ 3.6 Mpc [z|, we find 
a mass of 



Jiyf fOV 

JU DM,M87 



(0.75 ± 0.08) x lO lci M , (10) 



which agrees with the result of AFT (MJJ M Mg7 



f fov 



10 Mq) to within a factor of Vfov.cone/Vfov, 



prism 



tt/4. 



IV. THE DETECTABILITY OF v a DECAYS: 
GALAXIES VS. CLUSTERS 

In Fig. 2] we compare the diffuse spectrum of An- 
dromeda to the spectrum of Virgo A and determine the 
v s decay signals (in Counts/sec/keV) that would be pro- 
duced by each object for selected values of m s . To calcu- 
late the v s decay fluxes, we assume that sterile neutrinos 
comprise all of the dark matter, i.e., O s = f^DM = 0.24, 
and evaluate Eqn. © based on the distances and dark 
matter masses given in Table I for each object. Based 
on the total count rates and flux measurements of An- 
dromeda [64[ and Virgo A |j6(j, we divide by factors 

Of CM31 = g 3 x 10 -12 erg cm -2 Ct -1 and C M87 = 

7.0 x 10~ 12 erg cm~ 2 Ct _1 to convert the sterile neutrino 
fluxes (Eqn. to count rates. To realistically simulate 
the detected "line" fluxes in Fig. ^ we use a Gaussian 
centered at E 1 , s = m s /2 with a FWHM of AE = E 7jS /30 
(a conservative estimate of the energy resolution of the 
XMM EPIC detector 1 ). Doing so distributes ~ 72% of 
the signal over an energy range of ~ AE and converts 
the v s decay count rates to the same units as those of 
the measured spectra, Counts/sec/keV. 



1 http:/ /heasarc. nasa.gov/docs/xmm/xmm. html 
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TABLE I: Here we summarize the properties of Andromeda 

SUES] and Virgo A |H H3 |H [fl . In r0WS ^ " ^' We 
show, respectively, the distance to each object, the angular 
radius 9f ov of the XMM-Newton field of view (fov) of each 
observation, our estimates of the dark matter masses probed 
within each fov, the XMM exposure times (in kiloseconds: 
ks), and the (95% C.L.) upper bounds on m s . 



Galaxy Name 


Andromeda (M31) 


Virgo A (M87) 


Distance (Mpc) 


0.78 ± 0.02 


15.8 ±0.8 


#f ov (arcminutes) 


5.0' 


8.5' 


Mg£/10 U M S 


0.13 ±0.02 


75 ±8 


^cxp (ks) 


34.8 


25.9 


m a (keV) (95% C.L.) 


3.5 


8.2 



bound determined, through a similarly conservative pro- 
cess, from the spectrum of Virgo A (m s < 8.2 keV, 95% 
C.L.) 33] . The decay signatures corresponding to both 
of these limits are shown in Fig. ^ For comparison, we 
also show the decay signatures that would be present if 
Andromeda's dark matter halo were composed of sterile 
neutrinos with a mass of m. = 6.3 keV (the estimated 
Virgo A + Coma limit l36l Ir3l| . see Sec. I) or m s = 8.2 
keV (the Virgo A limit |33jp. The decay peak correspond- 
ing to 8.2 keV sterile neutrinos has a much more pro- 
nounced appearance in Andromeda than in Virgo A for 
two reasons: (1) the measured (background) spectrum 
of Andromeda is almost two orders of magnitude lower 
than that of Virgo A, yet, (2) regardless of the value of 
m s , is comparable to (Eqn. EJ: 



M31 



To determine the mass limit imposed by the diffuse 
spectrum of Andromeda, we evaluate the v s decay signal 
as described above for increasing values of m s until we 
reach the first statistically significant (> 2<7f) departure 
from the measured XMM spectrum. To be conservative, 
we use the largest fluctuations in the spectral data rela- 
tive to a smooth (power-law) fit to define the statistical 
significance (lcf) of the v s decay signal. The resulting 
2<7f limit is much more significant than the 95% C.L. 
defined by the formal statistical errors on the measured 
points. However, because features in the spectra may 
already reflect sterile neutrino decay and/or atomic line 
emission, detector backgrounds, etc. and are generally of 
uncertain origin (at least in the absence of detailed mod- 
eling), we argue that they are the appropriate gauge of 
statistical significance; i.e., for an upper bound on m s 
to be taken seriously, the corresponding decay signature 
should be large compared to any such features. In prac- 
tice, the limits we determine with this method roughly 
correspond to the lowest mass for which the decay sig- 
nal more than doubles the astrophysical background in 
a particular energy bin. In other words, to invalidate 
our upper bound the entire astrophysical background flux 
would have to vanish only in the bin in which we set our 
limit. 

In the case of Andromeda, the subtraction of compara- 
ble total and discrete point source emission within the 5' 
extraction region, L x ,tot = (2.2 ±0.2) x 10 39 erg s" 1 vs. 
L x ,pt = (2.0±0.1) x 10 39 erg s _1 , leads to a factor of two 
uncertainty in the normalization of the diffuse spectrum: 
£ x ,diffuse ^ (0.2 ± 0.2) x 10 39 erg s" 1 H. To account 
for this, we require the decay peak associated with the 
Andromeda mass limit to be at least four times the astro- 
physical background in the relevant energy bin. Based on 
this even more stringent criterion, i.e., 95% C.L.= 4<j[, 
the diffuse spectrum of Andromeda enables us to set an 
upper bound of m s < 3.5 keV (95% C.L.). Because of the 
rapidly increasing decay signal (oc m 3 ' 374 ) and rapidly 
falling background, marginally larger values of m s are 
excluded at even higher confidence levels. We note that 
our limit is significantly more restrictive than the upper 



cj)M87 



n 2 M iov 
^M87 -"^DM,M31 

r>2 II -f fov 
^M31 JW DM,M87 



0.71. 



(11) 



If the backgrounds of Andromeda and Virgo A were 
energy- independent, we would expect to be able to im- 
prove the Virgo A limit by a factor of ~ 100 3 ~ 4 
(Eqn. EJ), to m s iS 2 keV. However, because the X-ray 
spectrum of Andromeda rapidly increases at low ener- 
gies, as galactic and cluster X-ray spectra typically do, 
we do not quite reach this rough expectation. 



V. EXCLUSION LIMITS IN THE m a - sin 2 20 
PLANE 

To determine the region of the m s — sin 2 29 plane 
(Fig. |2J that is excluded by the X-ray spectra of An- 
dromeda and Virgo A, we must generalize the mass limits 
shown in Fig. ^ (which were derived under the assump- 
tion that fi s = 11dm = 0.24) to arbitrary values of il s . 
One approach is to derive a smooth exclusion region by 
scaling the total rate of v s radiative decays from a halo 
of mass A/dm, 



r s ,tot — 



Mi 



DM 



r, oc m 4 sin 2 20. 



(12) 



The key to this 'T s tot— scaling method", which is used 
in |23) LM > i s t° base the scaling on a v s decay peak 
(evaluated at fl s = 0.24) that is large compared to both 
the formal statistical errors on the flux as well as any 
fluctuations in the spectral data, as we have done. 

For each mass limit associated with such a peak, 
iVlim) we determine the mixing angle, sin 2#n s (via 
Eqn. |SJ, that corresponds to fl s = Odm = 0.24. By 
restricting alternative values of m s and sin 2 29 so that 
r s ,tot will remain fixed at r s , to t(w Si ii m , sin 2 20nj, we ar- 
rive at the following scaling of each mass limit through 
the m s — sin 2 26> plane: 



< 



sin 2 26> 
sin 2 26»o 



-1/4 



(13) 
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E [keV] 



FIG. 1: Here we compare the detectability of v s decays in 
Andromeda |64| and Virgo A |33l lfi(1|. The first statistically 
significant (4<Tf) v a decay peak relative to the measured spec- 
trum of Andromeda occurs at _E 7 , S = m s ,iim/2 = 1.75 keV, 
which excludes m a > 3.5 keV (95% C.L.). According to the 
analysis presented in [33|. the spectrum of Virgo A excludes 
m s > 8.2 keV (95% C.L.), which would produce a decay sig- 
nature like the dashed histogram. Because Andromeda would 
produce a similar v s decay signal to Virgo A (Eqn. [Till , but 
over a much smaller background, the prospective decay sig- 
nature of 8.2 keV sterile neutrinos in Andromeda is enormous 
by comparison. As an intermediate case, we also show what 
the decay peak associated with a 6.3 keV sterile neutrino, the 
estimated Virgo A + Coma mass limit (3ti loH , would look 
like in Andromeda. The vertical (la) error bars reflect the 
Poisson statistics of the signal and background count rates 
measured during each observation. 



In the case of Andromeda, we find 



rn , 



< 3.5 keV 



srrT20 



.94 x 10" 



1.91 keV 



sin 2 26> 

TrF 7 " 



-1/4 



-1/4 



(14) 



which excludes the region bounded by the dashed line 
(and above) in Fig. El 

There is an alternative and more accurate approach 
|6l| through which a limit on m s can be established more 
directly from the spectral data. For each value of m s , one 
determines the value of sin 2 26* for which <& X:S (m s , sin 20) 
(Eqn.^J exceeds the measured spectrum at E^_ s = m s /2 
by some threshold (4<7f in our case). We apply this "di- 
rect data method" to the (0.5—12 keV) diffuse spectrum 
of Andromeda. Above 5 keV, where the spectrum is dom- 
inated by instrumental background, we fix the amplitude 
at 0.003 Counts/sec/keV, which is well above the steeply 




sin 26 

FIG. 2: Here we present constraints on m s as a function of 
mixing angle, sin 2 2#, assuming that all dark matter is com- 
prised of sterile neutrinos. To facilitate comparisons, we adopt 
many of the conventions used by Abazajian and Koushiappas 
[3t^ . For L = 0, the thick, solid line corresponds to Q s = 0.24 
(Eqn.QjJ, while the shaded region to the right corresponds to 
Q, s > 0.24. Three density-production relationships associated 
with Q s = 0.3 and L 3> 10~ 10 are also shown [361 ]. The two 
previous direct radiative decay (v s —* ^ e . M , T + 7) u ppe r limits 
(both 95% C.L.) are based on measurements |55l l5fl I57I l5S| 
of the Cosmic X-ray Background [HIJ and XMM observa- 
tionsEglEl E3| of Virgo A (M87) and the Coma cluster 
[33ll36Ll6l| . The most stringent direct limits, from the present 
work (also 95% C.L.), are based on XMM observations of the 
Andromeda galaxy |64j . The region bounded by the dashed 
line is excluded by the T Sito t —scaling method" , while the re- 
gion above the solid, slightly jagged line is excluded by the 
more accurate "direct data method" (see Sec. 0. The in- 
direct lower limits (all 95% C.L.) labeled Lya(l) and Lyq(2) 
were derived in Ref. [33^ . while Lya(3) was derived in Ref. |33 |. 
Sterile neutrinos that occupy the horizontally hatched region 
could explain pulsar kicks |50l Isil Isl. I5H l54j . 



falling data |6J|. The conservative (95% C.L.) exclusion 
region that results is bounded by the slightly jagged line 
in Fig. An approximate fit to the boundary of this 
region is 



< 2.1 keV 



sin 2 26» 
10- 7 



-0.213 



(15) 



which is roughly parallel to the boundaries of the Cluster 
[361 Ifih] and CXB [5J| regions that were derived using 
essentially the same method (see Fig. |2 and Eqns. [HjI 
and IT7I below j. 

In addition to our Andromeda bounds, two previous 
radiative decay limits are also shown in Fig. EI Boyarsky 
et al. (5t| constrained m s using XMM-Newton (0.5—12 
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keV) [HIH and HEAO-I (3-60 keV) HHH measure- 
ments of the CXB. The resul ting limit (m s < 9.3 keV at 
fl s = 0.24) follows the trend |Hflj : 

m s < 4.1 keV j ^— j (16) 

According to the analyses in Refs. the XMM ob- 

servations of Virgo A (M87) [6(| and the Coma cluster 
[HIH demand: 

to, < 3.4 keV (^) . (17) 

Boyarsky et al. [7lJ recently determined a limit on 
sterile neutrino decays based on XMM observations of the 
Large Magellanic Cloud (LMC). We note that the signals 
shown in their spectra are very weak compared to back- 
grounds. Nevertheless, their claimed exclusion region is 
less restrictive than ours. Abazajian and Koushiappas 
[36| have also questioned the robustness of the LMC limit 
and all constraints based on dwarf/satellite galaxies due 
to the large uncertainties in their dark matter distribu- 
tions. 

The indirect constraints found in Ref. which 
we have also reproduced in Fig. |2 were derived using 
small scale clusterin g da ta from a variety of CMB mea- 
surements i[WMAP fijl Ell , CBI [13, Boomerang [U, 
ACBAR H3, VSA H2), the SDSS 3-D galaxy power 
spectrum, P g (k) p3| . the linear matter power spectrum 
inferred from Ly-a absorption in the SDSS quasar cata- 
log [2^, [2(| and from high-resolution observations of the 
Ly-a forest [271 Esl ]. When combined, the CMB measure- 
ments, SDSS 3-D P g {k), and SDSS Ly-a forest observa- 
tions set a lower bound (Lya(l)) of |3J| : 

m s >1.7keV (CMB + P g (k)+ Ly-a 95% C.L.). (18) 

The high-resolution Ly-a forest data could potentially 
provide an even stronger lower limit on to,, were it not 
for the (15% — 30%) systematic err ors |29j l30l l33| . In the 
limit of 15% (Gaussian) systematic I27ll28l uncertainties, 
combining the high-resolution data of 28j with the three 
data sets used to derive Eqn. ((TBI) yields (Lya(2)) |33j : 

m s > 3.0 keV (Previous + HR Ly-a 95% C.L.). (19) 

Recently, Seljak et al. j^ij recalculated this lower bound, 
finding m s > 14 keV (95% C.L.), which we also show in 
Fig- ID (Lya(3)). As discussed in the introduction, these 
results are still quite new and the subject of some con- 
troversy |33, |36| , but the basic agreement between these 
findings and the independent analysis of Viel et al. [3{| 
(to s > 10 keV 95% C.L.), strengthens the case for such 
a restrictive, indirect limit. In any case, it is valuable 
to improve both the direct and indirect constraints sep- 
arately. 



VI. CONCLUSIONS 

Warm Dark Matter (WDM) models of structure for- 
mation may more easily explain the low level of observed 
small scale clustering than standard ACDM galaxy for- 
mation scenarios. In this paper, we used the diffuse X-ray 
flux of the Andromeda galaxy (M31) [64[ to improve the 
radiative decay (y a — > v c ,^,,t + 7) upper limits on the 
mass m s of sterile neutrino WDM. In the context of the 
model described in 0j El H2> EH > our analysis of the 
diffuse spectrum of Andromeda requires m s < 3.5 keV 
(95% C.L.). Because of the rapidly increasing sterile 
neutrino decay signal (oc to;?' 374 ) and rapidly falling back- 
ground, larger values of m s are excluded at much higher 
confidence levels. As a case in point, we demonstrated 
that the decay signature associated with the upper limit 
(to s < 8.2 keV; 95% C.L. |33|) inferred from the spec- 
trum of Virgo A (M87) would be enormous relative to 
the significantly reduced astrophysical background of An- 
dromeda. When we combine our direct constraint with 
the most conservative, indirect lower limit set by mea- 
surements of small scale clustering in the CMB, the SDSS 
3-D galaxy power spectrum, and the Ly-a forest [3i| , we 
find that m s is restricted to the narrow range 

1.7 keV < to s < 3.5 keV (95% C.L.), (20) 

in the L ~ 10~ 10 ~ case. If the sterile neutrinos oc- 
cupy this tiny window, they could still be the dark matter 
and generate pulsar kicks jS^, EH HE El > but most 
of the parameter space remains viable only if the lepton 
asymmetry is very large: L ;§> 10~ 10 . Indeed the cor- 
roboration of the Seljak et al. (to s > 14 keV; 95% C.L.) 
and Viel et al. (m s > 10 keV; 95% C.L.) Lya constraints 
[3 EH strongly suggest that the standard L = produc- 
tion scenario of Abazajian et al. [I2I l3ll l32 | is no longer 
viable. 

A point source-subtracted XMM and/or Chandra spec- 
trum of Andromeda (and/or other nearby, massive yet 
quiescent spiral galaxies) extracted from within a radius 
of > 50' — 500' would probe > 10—100 times more dark 
matter than the observations we have considered here 
|64|. For a sufficiently large extraction region, the en- 
hancement in the predicted v s decay signal should be 
large enough compared to the additional diffuse, astro- 
physical background enclosed to provide sensitivity to 
sterile neutrino masses as low as m s = 1.7 keV. Detailed 
modeling of astrophysical emission and instrumental re- 
sponse, which we have not attempted in the present pa- 
per, might enable us to set even more stringent con- 
straints by allowing us to relax our limiting criterion 
without sacrificing robustness. In addition to provid- 
ing an even more restrictive bound on m s , this initia- 
tive would establish a much more comprehensive picture 
of the X-ray point source population and diffuse emis- 
sion of Andromeda (and/or other nearby spirals). The 
proposed observation(s) would therefore greatly benefit 
both the X-ray astronomy and astroparticle physics and 
cosmology communities. 
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At the same time, analyses of new and existing small 
scale structure data and the cutoff scale in the linear mat- 
ter power spectrum that they do (or do not) reveal will 
provide steadily improving lower bounds on m s . The con- 
vergence of these complementary constraints should lead 
to the detection or exclusion of sterile neutrino WDM in 
the near term. 
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